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摘 要 : 如 何 协调 区 域 社会 经 济 和 人 口 .资源 环境 发 展 ,已 成 为 生态 环境 建设 面临 的 紧迫 任务 ,识别 特定 地 区 的 脆 
弱 性 则 是 实现 有 效 生态 系统 管理 的 前 提 。 以 京 津 风沙 源 一 一 浑 善 达 克 沙 地 为 研究 区 ,采用 "暴露 -敏感 -适应 ” 
(VSD ) 模 型 建立 评价 指标 体系 ,结合 Moran 's 7 指数 与 主 成 分 分 析 , 对 2000 一 2019 年 该 沙 地 生态 脆弱 性 进行 评价 和 
驱动 机 制 分 析 。 结 果 表 明 : 浑 善 达 克 沙 地 生态 脆弱 性 20 a 间 呈现 增加 的 趋势 ;其 生态 脆弱 性 空间 上 呈 集 聚 现象 , 表 


现 为 西部 高 -高 聚集 .东部 低 - 低 聚集 的 分 布 格局 , 且 具 有 明显 的 空间 相关 性 ;生态 脆弱 性 的 时 空 分 布 的 异 质 性 是 由 


不 同 程度 的 人 类 活动 与 气候 因素 的 合 加 效应 造成 的 。 该 研究 为 浑 善 达 克 沙 地 的 生态 建设 和 生态 恢复 提供 理论 依 


据 , 并 提出 因地制宜 的 治理 建议 和 措施 。 
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在 全 球 气候 变化 及 人 类 活动 双重 作用 下 ,土地 
荒漠 化 .生物 多 样 性 丧失 .生态 系统 退化 等 一 系列 
生态 环境 问题 加 剧 "。 如 何 协调 区 域 社 会 经 济 和 人 
O ,资源 环境 发 展 ,已 成 为 生态 环境 建设 面临 的 紧 
迫 而 又 艰巨 的 任务 ”“。 科 学 认 知 评价 生态 环境 状 
态 是 进行 有 效 的 生态 管控 的 前 提 之 一 ,识别 特定 
地 区 的 生态 脆弱 性 则 是 应 对 这 些 变化 的 关键 ”。 

生态 脆弱 性 指 “ 随 着 时 间 和 空间 的 推移 ,生态 
系统 调节 其 对 压力 源 的 反应 的 潜力 ”。 作 为 可 持 
续 科学 的 一 部 分 ,国内 外 学 者 在 理论 和 方法 等 方面 
对 生态 脆弱 性 开展 了 大 量 工作 ,评价 区 域 涉及 城市 
和 工矿 开发 “自然 保护 区 ”黄土 高 原 ”\ 生 态 脆 
35 pc ”等 地 区 , 且 对 我 国 北方 草原 区 沙 地 生态 系 
统 的 研究 相对 较 少 。 在 评价 方法 上 ,通常 从 上 自然 和 
人 类 活动 两 方面 选取 多 个 指标 组 成 评价 体系 ”, 考 
虑 到 评价 体系 中 主观 指标 可 能 会 导致 评 佑 结果 的 
不 确定 性 等 问题 ,伴随 着 数据 的 获取 手段 的 增加 ， 
学 者 们 逐渐 倾向 于 使 用 由 遥感 手段 获取 的 客观 指 
标 "”。 为 满足 对 评价 指标 和 信息 数据 的 高 效 、 系 统 
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的 组 织 需 求 "“ ,充分 体现 “人 一 地 ”耦合 关系 ,生态 
脆弱 性 理论 体系 逐渐 从 以 “因素 一 结果 ”的 方法 体 
系 ,如 “压力 一 状态 一 响应 ”模型 (Pressure State Re- 
sponse ,PSR)”、“ 驱 动力 一 压力 一 状态 一 影响 一 响 
应 一 管理 ”模型 (Driving-force Pressure State Impact 
Response Management, DPSIRM) "等 ,向 强调 “ 暴 
FE 敏感 性 和 适应 性 "的 体系 转变 “” ,如 “暴露 一 敏 
感性 一 应 对 能 力 ”(Exposure Susceptibility Coping Ca- 
pacity ,ESC) 模 型 中“ 敏感 性 一 恢复 力 一 压力 度 ” 模 
型 (Sensitivity Recovery Pressure, SRP)! 等 。 其 中 ， 
Polsky 等 "在 后 者 体系 上 发 展 了 VSD 模 型 (Vulnera- 
bility Scoping Diagram) 评 价 流程 , 因 评 价 指标 建立 过 
程 逐 级 递 进 流程 规范 清晰 ,得 到 较 广泛 应 用 ?”。 
近年 来 , 随 着 数据 获取 手段 .研究 方法 和 研究 
理论 的 不 断 发 展 ,对 生态 脆弱 性 时 空 动态 性 和 驱动 
机 制 逐渐 进入 研究 者 的 视野 ,应 用 半 方 差 函 数 、 空 间 
自 相 关 等 研究 生态 脆弱 性 的 空间 异 质 性 ”或 识别 
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变 及 动因 ,应 用 地 理 探测 器 分 析 陕 西 省 和 江西 省 万 
安县 生态 脆弱 性 时 空 演 变 及 驱动 力 " 56, PATE, 
随 着 研究 时 间 跨 度 的 增加 ,导致 生态 脆弱 性 的 驱动 
力 可 能 会 发 生变 化 ,如 气候 变化 、 环 境 保 护 政策 的 
实施 .土地 利用 方式 的 改变 等 。 对 于 导致 生态 脆弱 
性 变化 的 自然 因素 和 人 类 活动 因素 的 驱动 强度 是 
否 发 生变 化 ,目前 研究 较 少 。 

因 京 津 风沙 源 等 生态 修复 工程 的 实施 ,2000 年 
后 浑 善 达 克 沙 地 的 生态 环境 有 所 改善 WW。 但 因 草 
AFE 、. 过 度 放牧 等 不 合理 的 土地 利用 与 气候 暖 干 
化 ,其 脆弱 的 生态 系统 更 易 受到 干扰 ,甚至 在 某 些 
地 区 沙漠 化 再 度 发 展品, 导致 该 地 区 可 持续 发 展 前 
景 不 容 乐 观 2 32。 因此 ,本 研究 以 浑 善 达 克 沙 地 为 
人 研究 区 ,基于 VSD 模 型 构建 沙 地 生态 脆弱 性 评 佑 指 
标 体 系 ,通过 评估 该 地 区 2000—2019 4F 5 HAE AS Jf 
弱 性 ,揭示 其 生态 脆弱 性 的 时 空 分 布 特征 ,探讨 导 
致 其 脆弱 性 变化 的 驱动 机 制 , 以 期 为 该 范围 内 的 生 
态 建设 和 生态 恢复 提供 理论 依据 ,提出 因地制宜 的 
治理 措施 和 建议 。 


1 数据 和 方法 


1.1 研究 区 概况 
浑 善 达 克 沙 地 位 于 内 蒙古 中 部 锡林郭勒 草原 
南端 ,是 我 国 著名 的 四 大 沙 地 之 一 ,也 是 困扰 京 津 
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责 的 沙 侍 的 源头 之 一 ,该 区 生态 系统 极其 脆弱 。 坐 
标 范围 为 112.0°~117.5°E、42.0°~43.5°N 之 间 , 东 西 
长 450 km ,平均 海拔 1200 m; 属 中 温带 大 陆 性 气候 ， 
年 均 温 1.5 "C ;年 降雨 量 约 为 360 mm ,全 年 降雨 量 的 
80% 以 上 集中 在 7 一 9 月 。 
1.2 生态 脆弱 性 指标 体系 构建 

参考 相关 研究 的 基础 上 2 ,选取 13 个 指标 构 
建 浑 善 达 克 沙 地 生态 脆弱 性 评价 体系 ( 表 1) ,分 别 
对 2000 年 .2005 年 .2010 年 .2015 年 .2019 年 生态 脆 
弱 性 进行 分 析 。 

土地 利用 / 窗 被 数据 ,基于 中 国 科 学 院 地 理科 学 
与 资源 研究 所 1:100000 土 地 利用 数据 ,构建 了 2000 
年 .2005 年 .2010 年 .2015 年 浑 善 达 克 沙 地 土地 利用 
数据 集 ,其 中 2019 年 数据 来 自 对 Landsat 8 OLDER 
数据 的 土地 利用 更 新 解 译 。 

气象 数据 来 源 于 “中 国 地 面 气候 资料 日 值 数 据 
集 (V3.0)”(http:/data.cma.cn/)。 利 用 浑 善 达 克 沙 地 
周边 气象 站 点 ,获得 年 均 降 水 、 生 长 季 平 均 降 水 间 
隔 、 极 端 最 高 气温 天 数 (>33 C ) 等 指标 ,在 ArcGIS 
10.3 中 采用 克 里 金 插值 法 生成 棚 格 图 层 。 而 年 降水 
变化 率 , 则 选择 研究 时 有 段 前 后 共 3 a, 以 Theil-Sen 
median 趋势 表示 (公式 见 文献 [33])。 

地 温 采 用 MODIS MOD11A2 地 表 温 度 产 品 
(2000—2019 年 ,时 间 分 辩 率 8 d 空间 分 辩 率 250 


R1 生态 脆弱 性 评价 指标 体系 及 权重 


Tab. 1 The ecological vulnerability evaluation indices system and weight 


标 层 (权重 ) 要 素 层 (权重 ) 指标 层 * 指标 方向 指标 权重 

暴露 度 BI 人 口 社 会 压力 A OREA km?) + 0.0384 
(0.4934) (0.2470) ASS SHOT BY em? A) - 0.0488 
TERR GE km?) 十 0.1598 
B2 干旱 灾害 地 温 /C 十 0.0156 
(0.2465) 年 降水 变化 率 /(mm':an) + 0.1027 
生长 季 降 水 间隔 /d + 0.1112 
极端 最 高 温 天 数 /d + 0.0170 
敏感 性 B3 植被 覆盖 NDVI - 0.0940 
(0.3108) (0.1726) 景观 多 样 性 B 0.0245 
景观 破碎 度 十 0.0541 
B4 气 候 年 均 降水 量 /mm - 0.1382 

(0.1382) 
适应 能 力 B5 自我 调节 NPP/(g-m?+a") 一 0.1305 
(0.1958) (0.1958) NDVI 变 化 率 " 0.0653 


内 , 即 值 越 大 对 生态 脆弱 性 影响 就 越 大 。 


注 :a 表示 所 有 棚 格 数据 空间 分 辩 率 为 250 m, 横 轴 墨 卡 托 投 影 (WGS 1984 UTM Zone 49N) ,空间 分 辨 率 250 m。“+” 表 示 该 指标 越 大 越 增加 系统 
脆弱 ;…- "表示 该 指标 越 大 有 利于 减少 系统 脆弱 ;根据 指标 对 脆弱 性 影响 的 情况 ,分别 对 正 ,. 逆 指标 进行 标准 化 ,将 指标 值 归 一 化 到 0~1 范 十 
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m, https://ladsweb. modaps.eosdis.nasa.gov/search/ ) ,在 Y Yu, A NE. 
ArcGIS 10.3 软 件 中 采用 5 一 9 月 生长 季 平 均值 合成 。 pes (2) 


归 一 化 植被 指数 (NDVI) ,采用 MODIS 
MOD13Q1(2000 一 2019 年 ,时 间 分 辩 率 16 d 空间 分 
3953€ 250 m) ,在 ArcGIS 10.3 软件 中 采用 最 大 值 合成 
而 得 。NDVI 变 化 率 ,以 3 a 为 时 间 段 ,采用 Theil- 
Sen median 趋势 分 析 计算 。 

植被 净 初 级 生产 力 (NPP) ,采用 CASA 模 型 9 ， 
该 模型 考虑 太阳 辐射 .温度 、 水 分 等 胁迫 因子 对 植 
被 NPP 的 影响 ,计算 公式 见 文献 [35]。 

景观 多 样 性 和 景观 破碎 度 ,以 栅 格 化 后 的 5 期 
土地 利用 数据 为 对 象 ,在 Fragstats 4.2 软件 中 计算 
获得 。 

AD Bt .牲畜 数量 等 社会 经 济 统计 数据 来 源 
于 历年 《锡林郭勒 盟 统计 年 鉴 >。 人 口 密 度 为 各 乡 
镇 (或 苏 木 ) 的 总 人 口 /总 面积 ;人 均 草 地 面积 为 各 乡 
镇 (或 苏 木 ) 的 草地 总 面积 /总 人 口 ;牲畜 数量 图 层 则 
基于 NDVI 加 权 的 面 插值 方法 获得 *。 

指标 权重 采用 层次 分 析 法 (AHP) 计 算 获 得 ,应 
用 yaahp10.3 软件 (https://www.metadecsn.com/yaahp/)， 
采用 1~9 标 度 法 分 别 对 目标 层 和 各 要 素 层 的 指标 进 
行 重要 性 两 两 比较 打分 ,构建 判断 矩阵 ,并 通过 一 
致 性 检验 (检验 值 分 别 为 0.052、0.100、0.091 和 
0.000) ,表明 各 要 素 权重 结果 可 信和 度 较 高 ( 表 1)。 

13 生态 脆弱 性 计算 
生态 脆弱 性 指数 (EVI) 公 式 如 下 : 


EVI =X w,X, (1) 
i=l 


式 中 :EVI 为 生态 脆弱 性 指数 ; w, 为 第 ;个 指标 权 
E; X, 为 第 ;个 指标 图 层 。 取 5 期 评价 结果 的 自然 
断 点 法 分 级 (Natural breaks ) 的 均值 ,作为 最 终 的 分 
级 结果 2 : 微 度 脆弱 (EVI<0.34) 、 轻 度 脆弱 (0.34 < 
EVI« 0.39) .中 度 脆弱 (0.39 < EVI< 0.44) .重度 脆弱 
(0.44 < EVI< 0.50) ,极度 脆弱 (EVI > 0.50). 
14 空间 自 相关 分 析 

为 精确 表达 浑 善 达 克 沙 地 生态 脆弱 性 空间 特 
征 ,在 研究 区 布设 2.5 kmx2.5 km 网 格 ,每 个 网 格 即 
生态 脆弱 性 单元 区 , 取 值 为 该 网 格 所 包围 区 域内 胸 
弱 性 的 均值 。 在 ArcGIS 10.3 软件 中 ,选用 全 局 和 局 
部 Moran's7 指 数 ,分 别 对 该 沙 地 5 期 生态 脆弱 性 指 
数 的 空间 差异 及 集聚 特征 进行 分 析 , 公 式 如 下 : 


式 中 : w 为 网 格 i 与 网 格 j 之 间 的 空间 权重 矩阵 ; 
x, x 为 网 格 的 属性 值 ; n 为 空间 单元 数量 ; S 为 方 
差 。Moran”s 1 值 的 区 域 为 -1~1。Moran's />0 时 , 表 
示 空 间 正 相关 性 ,其 值 越 大 ,空间 相关 性 越 明 显 ， 
Moran’ s 7<0 表 示 空 间 负 相 关 性 ,其 值 越 小 ,空间 差 
异 越 大 ,否则 ,Moran's /=0, 空 间 呈 随机 性 。 
1.5 生态 脆弱 性 驱动 因素 变化 分 析 

为 了 分 析 人 类 活动 对 该 地 区 生态 脆弱 性 的 影 
响 , 应 用 5 期 土地 利用 图 ,以 各 土地 利用 类 型 为 统计 
单元 ,使 用 ArcGIS 10.3 的 Zonal statistics as Table 工 
具 , 获 得 不 同 研究 时 段 各 土地 利用 类 型 的 生态 脆弱 
性 均值 ,以 表征 人 类 活动 的 影响 。 

采用 主 成 分 分 析 法 (PCA) ,对 5 期 该 沙 地 生态 
脆弱 性 评价 体系 的 13 个 指标 进行 分 析 , 识 别 影响 其 
生态 脆弱 性 的 驱动 因素 及 其 变化 特征 。 其 中 ,驱动 
因素 的 强 弱 以 载荷 系数 表示 ,该 系数 类 似 于 回归 系 
数 , 值 越 大 表示 该 驱动 因子 对 生态 脆弱 性 影响 越 
大 。 主 成 分 分 析 在 SPSS 22.0 软件 中 完成 。 


2 结果 与 分 析 


21 生态 脆弱 性 时 空 变化 分 析 

从 时 间 变 化 看 ,2000 一 2019 年 浑 善 达 克 沙 地 生 
态 脆弱 性 由 中 度 脆 弱 过 渡 到 重度 脆弱 ,总 体 脆弱 性 
旺 增加 的 趋势 ,平均 脆弱 性 指数 分 别 为 :0.407、 
0.411 ,0.469 .0.462 .0.497。 从 空间 分 布 来 看 ,该 沙 
地 脆弱 性 呈现 西部 区 域 最 高 .中 部 区 域 次 之 、 东 部 
区 域 最 低 的 分 布 格局 (图 1)。 极 度 和 重度 脆弱 区 分 
布 在 西部 区 域 , 从 2000 年 开始 ,两 类 脆弱 区 逐渐 向 
东部 扩张 ,极度 脆弱 区 由 2000 年 的 2.0% 扩 展 到 
2019 年 的 47.5% ,重度 脆弱 性 分 布 区 域 变 化 与 之 类 
似 ,面积 在 2010 年 达到 最 大 值 (70.9%)。 除 2015 年 
外 ,20 a 间 该 沙 地 中 度 脆弱 区 面积 呈现 减少 趋势 , 面 
积 比 由 48.3% 锐 减 到 16.6%。 轻 度 脆弱 区 由 该 地 区 
的 中 东部 向 东部 退缩 , 自 2000 年 的 35.7% 锐 减 到 
2019 年 的 2.4%。 
2.2 生态 脆弱 性 空间 自 相关 分 析 

2000—2019 年 浑 善 达 克 沙 地 生态 脆弱 性 具有 
明显 的 空间 自 相 关 , 全 局 Moran’ s 了 指数 分 别 为 
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注 :(a)~(e) 为 各 年 份 生态 脆弱 性 分 级 , (f) 为 各 年 代 各 级 生态 脆弱 性 面积 比 。 
1 2000 一 2019 年 浑 善 达 克 沙 地 生态 脆弱 性 空间 分 布 
Fig. 1 The spatial distribution of ecological vulnerability in Otindag Sandy Land from 2000 to 2019 


0.9833 , 0.9882 .0.9743 .0.9904 和 0.9887。 局 部 Mo- 
ran"s7 了 指数 呈 区 域 分 布 (图 2) , 即 高 -高 聚集 区 集中 
在 该 地 西部 , 低 - 低 聚集 区 则 集中 在 东部 区 域 。 
2000 年 ,高 -高 聚集 区 面积 占 比 为 10.2% , 主要 分 布 
在 苏 尼 特 右 旗 中 部 , 低 - 低 聚集 区 集中 在 浑 善 达 克 
沙 地 的 中 东部 区 域 ,面积 占 比 为 5.3%。2005 年 ,高 - 
高 聚集 区 向 该 沙 地 西部 收缩 ,面积 比 减少 到 7.8% ; 
低 - 低 聚集 区 向 东 收 缩 ( 面 积 比 43.7% )。2010 年 ， 
高 -高 聚集 区 面积 占 比 分 别 比 2005 年 增加 7.8%; 
低 - 低 聚集 区 面积 较 2005 年 减少 20.3% , 且 由 连续 
分 布 状态 转变 成 破碎 化 分 布 状态 (面积 比 43.7%)， 
表明 该 时 间 段 浑 善 达 克 沙 地 脆弱 性 增加 。2015 年 
与 2019 年 ,高 -高 聚集 区 扩展 的 趋势 得 到 缓解 ,面积 
占 比 分 别 为 11.7% 与 10.8% , 低 - 低 聚集 区 域 面积 
所 增加 ,面积 比分 别 为 40.7% 与 37.8%。 
2.3 不 同 土地 利用 类 型 的 生态 脆弱 性 特征 

对 2000—2019 年 不 同 土地 利用 生态 脆弱 性 排 
序 为 ( 表 2) :水 域 > 未 利用 土地 > 草地 > 居民 工矿 用 


地 > 耕地 > 林地 。 耕 地 、 林 地 与 居民 工矿 用 地 的 生态 
脆弱 性 20 a 间 变化 趋势 基本 相似 。 三 者 在 2005 年 
出 现 最 低 值 ,在 2010 年 出 现 最 高 值 ;2010 年 后 3 种 
土地 利用 类 型 的 脆弱 性 表现 出 下 降 或 持平 的 趋 
势 。 草 地 、 水域 和 未 利用 土地 的 生态 脆弱 性 在 
2000 一 2019 年 间 呈 增加 的 趋势 ,最 低 值 出 现在 2000 
年 (分 别 为 0.407 .0.419 .0.410) ,最 高 值 出 现在 2019 
年 (分 别 为 0.450 .0.466 .0.455 ) 。 
24 生态 脆弱 性 驱动 因素 分 析 

采用 主 成 分 分 析 识 别 影响 浑 善 达 克 沙 地 生态 
脆弱 性 的 驱动 因素 及 其 变化 特征 。 根 据 共同 因素 
可 以 解释 指标 变量 50% 变 异 量 的 原则 ,选取 因子 载 
Ti 5 0.71 ,此 时 因子 载荷 量 的 状况 "其 为 理想 ”7 
( 表 3)。 对 该 地 区 各 研究 时 段 第 一 主 成 分 可 知 , 生 
态 脆弱 性 主要 与 反映 人 类 活动 的 指标 (人 口 密度 、 
牲畜 数量 等 ) 气候 变化 指标 (年 均 降 水 .生长 季 降 
水 间隔 .年 降水 变化 率 和 地 温 等 ) .植被 状态 (NDVI 
和 NPP 等 ) 指 标 有 关 。 第 二 主 成 分 与 景观 多 样 性 有 
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注 :(a)~(e) 为 各 年 份 生态 脆弱 性 局 部 Moran s 7 分布 ,(f) 为 
自身 生态 脆弱 性 和 周边 水 平 都 比较 高 , 低 - 低 聚集 
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局 部 Moran’ s 17 不 同 聚 类 类 型 面积 比 ;高 -高 聚集 区 表示 
区 表示 自身 生态 脆弱 性 和 周边 水 平 都 比较 低 。 


图 2 2000 一 2019 年 浑 善 达 克 沙 地 生态 脆弱 性 局 部 Moran ' s 7 分 布 
Fig. 2 Moran’s / cluster diagram of ecological vulnerability of Otindag Sandy Land from 2000 to 2019 


表 2 2000—2019 年 不 同 土地 利用 类 型 的 生态 环境 脆弱 指数 
Tab. 2 Ecological vulnerability index of different land use types from 2000 to 2019 


土地 利用 类 型 (2019 年 面积 比 ) 2000 年 2005 年 2010 年 2015 年 2019 年 平均 值 
耕地 (3.18% ) 0.377 0.356 0.447 0.424 0.425 0.406 
林地 (0.83% ) 0.375 0.363 0.430 0.413 0.419 0.400 
草地 (77.96% ) 0.407 0.412 0.469 0.463 0.498 0.450 
水 域 (1.69%) 0.419 0.437 0.482 0.470 0.521 0.466 
居民 工矿 用 地 (0.31%) 0.386 0.371 0.455 0.435 0.446 0.419 
未 利用 土地 (16.03%) 0.410 0.417 0.472 0.469 0.508 0.455 
关 , 反 映 了 该 地 区 土地 利用 变化 情况 。 然 而 ,在 20 a 
3 讨论 


的 时 间 跨 度 中 ,不 同年 代 的 生态 脆弱 性 的 驱动 因素 
略 有 不 同 , 人 口 密度 .牲畜 数量 分 别 在 2010 年 .2015 
年 对 生态 脆弱 性 的 影响 减弱 ,人 均 草 地 面积 指标 也 
有 这 种 趋势 ;降水 因素 .NDVI 与 NPP 在 整个 研究 时 
段 内 均 与 生态 脆弱 性 关联 较 大 ,降水 的 多 少 .植被 
状态 的 好 坏 影 响 着 当地 的 生态 脆弱 性 。 景 观 多 样 
性 指标 在 2000 年 .2005 年 与 该 地 区 的 生态 脆弱 性 关 
系 较 大 。 


3.1 生态 脆弱 性 的 自然 和 人 类 活动 驱动 因素 变化 
浑 善 达 克 沙 地 生态 脆弱 性 受到 人 类 活动 与 气 
候 因 素 的 双重 影响 ,导致 生态 脆弱 性 时 空 分 布 异 
质 性 。 
(1) 从 不 同 土地 利用 类 型 的 生态 脆弱 性 变化 来 
看 ( 表 2) ,生态 脆弱 性 与 人 类 土地 利用 情况 有 关 ” 。 
水 域 和 未 利用 土地 表现 出 较 高 的 生态 脆弱 性 ,可 能 
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表 3 2000 一 2019 年 各 研究 时 段 生 态 脆弱 性 主 成 分 载荷 矩阵 
Tab.3 The load matrix of principal components from 2000 to 2019 


2 2000 4E. 2005 年 20104F. 2015 年 2019 年 
评价 指标 
PC1 PC2 PC1 PC2 PC1 PC2 PC1 PC2 PC1 PC2 
人 口 密度 0.76 0.77 
人 均 草 地 面积 0.86 0.71 
牲畜 数量 0.84 0.90 0.81 
也 温 0.79 0.77 0.81 0.80 
年 降水 变化 率 0.90 0.91 
生长 季 降 水 间隔 0.92 0.93 0.89 0.93 0.94 
极端 最 高 温 天 数 0.84 0.92 0.87 0.82 
NDVI 0.83 0.84 0.85 0.87 0.85 
景观 多 样 性 -0.71 -0.79 
景观 破碎 度 
年 均 降 水 0.96 0.96 0.96 0.95 0.95 
NPP 0.88 0.85 0.87 0.83 
NDVI 变 化 率 
贡献 率 /% 51.62 10.46 59.33 9.76 54.14 10.94 53.73 10.74 53.86 11.01 


注 : 表 中 仅 保留 了 | 因子 载荷 (factor loadings)l>0.71 的 因子 载荷 ; 正 负 向 指标 均 通 过 归 一 化 ,调整 为 指标 值 越 大 对 生态 脆弱 性 影响 越 大 。 


因 年 内 /年 际 间 降 水 和 蒸发 导致 水 域 面积 发 生变 化 ， 
裸露 的 滩地 出 现 较 高 的 植被 覆盖 度 及 未 利用 土地 
会 出 现 较 多 的 一 二年生 植 物 ,因而 表现 出 较 高 的 
生态 脆弱 性 ;林地 拥有 较 高 的 生态 系统 稳定 性 (或 
较 高 的 抗 外 界 干扰 能 力 ), 但 同时 20a 间 林地 的 生态 
脆弱 指数 从 2000 年 的 0.375 增 加 到 2019 年 的 0.419， 
造成 的 原因 多 为 榆树 疏 林 是 牧民 放牧 的 夏 营 地 , 放 
牧 导 致 榆树 幼苗 被 采 食 严重 ,榆树 种 群 更 新 受阻 、 
多 年 生 草 本 比例 减少 、 裸 沙 面 积 增加 ;2000 一 
2019 年 ,草地 发 生 退 化 ,高 .中 、 低 覆盖 草地 面积 4 
别 变 化 -422.82 km? -220.34 km, 1699.76 km ,草地 
的 生态 脆弱 性 较 高 , 且 在 2019 年 脆弱 性 指数 逐渐 增 
加 至 重度 脆弱 (0.498) ;耕地 因 人 工 灌 溉 等 因素 , 表 
现 出 比 草地 更 低 的 脆弱 性 。 

(2) 浑 善 达 克 沙 地 脆弱 性 的 空间 分 布 特征 亦 与 
气候 指标 有 关 。 从 东部 的 典型 草原 带 , 过 渡 到 西部 
的 荒漠 草原 带 , 降 水 量 逐 渐 减 少 , 导 致 生态 脆弱 性 
呈 聚 集 状 分 布 , 即 轻 度 脆弱 区 ( 低 - 低 聚集 ) 分 布 在 
该 地 区 的 东部 .重度 脆弱 区 域 ( 高 -高 聚集 ) 分 布 在 
西部 (图 2)。 加 之 我 国 北方 气候 暖 干 化 的 发 展 趋 
35 ,进一步 加 剧 了 该 地 区 生态 系统 脆弱 性 。 而 在 
临近 的 坝 上 地 区 , 因 降 水 作用 的 影响 其 生态 脆弱 性 
从 东 向 西 增强 中 ,以 及 在 黄土 高 原 神 东 矿区 生态 脆 
弱 性 也 受降 水 因素 的 主导 "* ,这些 地 区 生态 脆弱 性 


均 表现 出 与 本 研究 区 域 相似 的 变化 特征 。 

(3) 导致 该 地 区 生态 脆弱 性 变化 的 自然 和 人 类 
活动 因素 在 不 同 研 究 时 段 的 到 加 效应 存在 差异 。 
由 主 成 分 分 析 可 知 ( 表 3) ,反映 人 类 活动 的 人 口 密 
BE .牲畜 数量 和 人 均 草 地 面积 等 指标 在 2010 年 或 
2015 年 后 对 生态 脆弱 性 的 影响 减弱 。 随 着 生态 恢 
复工 程 和 生态 政策 的 实施 等 ”, 人 类 活动 增强 的 趋 
势 受 到 抑制 ,而 逐渐 凸显 气候 因素 对 该 地 区 生态 脆 
弱 性 驱动 的 主导 作用 。 对 2000 一 2015 年 锡林郭勒 
植被 NDVI 变 化 驱动 力 定 量 分 析 也 发 现 了 类 似 的 现 
象 ,NDVI 显 车 减少 区 域 ,气候 因素 的 驱动 作用 超过 
人 类 活动 因素 ,而 在 NDVI 显 著 增 加 区 域 ,两 者 关系 
则 相反 2 。 

3.2 政策 建议 

从 区 域 尺 度 考 虑 ,以 未 来 空间 管制 角度 出 发 ， 
浑 善 达 克 沙 地 东部 区 域 生态 脆弱 性 以 轻 度 .中 度 为 
主 , 生 态 管理 措施 应 侧重 于 草地 资源 合理 利用 与 生 
态 保护 并 举 ,避免 生态 脆弱 性 提高 ,充分 发 挥 " 京 津 
囊 " 生 态 服 务 供应 地 的 作用 ;而 该 沙 地 中 西部 地 区 
脆弱 性 则 以 重度 及 以 上 为 主 , 需 重 点 强化 生态 恢复 
工程 .生态 移民 减少 放牧 强度 等 ,恢复 生态 。 

从 局 地 尺度 人 类 活动 方面 来 说 ,过 度 放牧 导致 
该 沙 地 生态 系统 退化 ,进而 导致 生 态 脆弱 性 提高 ， 
如 何 减少 放牧 压力 是 解决 当地 生态 脆弱 性 持续 恶 
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Spatio-temporal patterns variation of ecological vulnerability in Otindag Sandy 
Land based on a vulnerability scoping diagram 


CHEN Zhenqi', ZHANG Jing", ZHANG Yilong, LIU Rui’ 


(1. College of Environment and Resources, Dalian Minzu University, Dalian 116600, Liaoning, China; 


2. School of Geography, Liaoning Normal University, Dalian 116029, Liaoning, China) 


Abstract: The question of how to coordinate social and economic development, population growth, resource 
consumption, and environmental protection has become an urgent topic for regional sustainability development. 
Identifying the ecological vulnerability (EV) of specific areas is the precondition to achieving effective ecosystem 
management for sustainability. However, the driving forces of EV, including climate change, environmental 
protection policies, and land use, may change over a long period of time, affecting the spatio-temporal patterns of 
EV. Taking Otindag Sandy Land as the research area, we used the “Vulnerability Scoping Diagram” framework 
to establish the evaluation index system. Using Moran’ s / index and principal component analysis (PCA), we 
evaluated the spatio-temporal distribution patterns of EV of this sandy land, and we also analyzed its driving 
mechanism from 2000 to 2019. The findings revealed that: (1) the EV of Otindag Sandy Land has exhibited an 
increasing trend over the past two decades; (2) the EV of this area presents an agglomeration distribution pattern 
with obvious spatial auto-correlation on spatial extent, in which the high-high cluster of EV is distributed in the 
western region, and the low-low cluster of EV is located in the eastern region; (3) the heterogeneity of the 
temporal and spatial distribution of EV is due to the synergistic effect of human activities and climate factors at 
different levels. According to the factor loadings matrix of PCA, after 2010 or 2015, the influence of indicators 
such as the population density, the number of livestock, and the per capita grassland area on the EV of this region 
has either weakened or been restrained by the implementation of ecological restoration projects and ecological 
policies. Meanwhile climate factors have gradually highlighted the leading role in driving EV in this region. This 
study provides a theoretical basis for the ecological construction and ecological restoration of Otindag Sandy 
Land, and puts forward some suggestions and measures for ecosystem management. At the regional level, from 
the perspective of future spatial planning, the EV in the eastern part of Otindag Sandy Land is mild to moderate, 
and ecological management measures should focus on harmonizing rational utilization of grassland resources 
with ecological protection. Meanwhile, the EV in the central and western regions of this area is severe to extreme, 
and it is thus necessary to strengthen ecological restoration projects and ecological migration to reduce grazing 
pressure on sandy land. At the local level, while comprehensively considering the demographic factors, pasture 
size, and grassland quality of herdsmen households, we should optimize grassland property rights and coordinate 
the interests of herdsmen with the goal of ecological conservation to achieve the fine administrative management 
“one household, one eco-management measures”. 


Keywords: ecological vulnerability; spatio-temporal pattern; driving forces; Otindag Sandy Land 


